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In this paper, we report results for the wave packet dynamics in a class of quasiperiodic chains 
consisting of two types of weakly coupled clusters. The dynamics are studied by means of the return 
probability and the mean square displacement. The wave packets show anomalous diffusion in a 
stepwise process of fast expansion followed by time intervals of confined wave packet width. Applying 
perturbation theory, where the coupling parameter v is treated as perturbation, the properties of the 
eigenstates of the system are investigated and related to the structure of the chains. The results show 
the appearance of non- localized states only in sufficiently high orders of the perturbation expansions. 
Further, we compare these results to the exact solutions obtained by numerical diagonalization. This 
shows that eigenstates spread across the entire chain for v > 0, while in the limit v — » ergodicity 
is broken and eigenstates only spread across clusters of the same type, in contradistinction to 
trivial localization for v = 0. Caused by this ergodicity breaking, the wave packet dynamics change 
significantly in the presence of an impurity offering the possibility to control its long-term dynamics. 

PACS numbers: 71.23.Ft, 71.15.-m, 72.15.-v 



I. INTRODUCTION 

Understanding the relations between the spectral prop- 
erties of a given Hamiltonian and the dynamics of wave 
packets which are governed by it, remains one of the ele- 
mentary questions of quantum mechanics that still poses 
significant challenges, further emphasized by the discov- 
ery of quasiervstals. While spectra of many Hamilto- 
nians decompose into a point-like part and an absolutely 
continuous part accompanied by bounded (localized) and 
unbounded (delocalized or extended) eigenstates, there 
exists a large variety of Hamiltonians whose spectra, for 
certain values of the parameters, are neither pure point- 
like nor absolutely continuous, nor a combination of both. 
In this case, the spectrum contains a singular continuous 
part, and eigenstates arc often found to be multifractal. 
Examples are Harper's model of an electron in the mag- 
netic fieldj^ the kicked rotator— £ as well as the Anderson 
model of an electron in a disordered medium^ In the case 
of an electron in a one-dimensional quasiperiodic system, 
as studied in this paper, many examples lead to spectra 
which are purely singular continuous as welli^jSjiS 

To address the above mentioned challenge, we investi- 
gate wave packet dynamics in one-dimensional quasiperi- 
odic chains by numerical simulation as well as pertur- 
bation theory, and relate the results to the hierarchical 
properties of these chains. The content of this paper is 
organized as follows: At first we introduce the construc- 
tion rule and structure of the chains in Sec.Ql] Sect ion Hill 
then focuses on numerical results for the time evolution 
of wave packets. To obtain a better understanding of the 
properties of wave packet spreading and localization, we 
apply perturbation theory in Sec. IIVI and further study 
the influence of an impurity on the wave packet dynamics 
in Sec. |V1 followed by a brief summary of our results. 



II. QUASIPERIODIC CHAINS WITH GOLDEN, 
SILVER OR BRONZE MEANS 



In this paper we study one-dimensional quasiperiodic 
systems constructed by the inflation rule 
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iterated a times starting from the symbol w, where the 
letter w denotes a weak bond and s a strong bond. We re- 
fer to the resulting sequence after a iterations as the ath 
order approximant C a with the length f a given by the 
recursive equation /„ = / a _ 2 + nf a -\ and / = fx = 1. 
Depending on the parameter n the inflation rule gener- 
ates different so-called metallic means, i.e. the lengths of 
two successive sequences satisfy the relation 
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where A is an irrational number with the continued frac- 
tion representation [n] — [n, n, n, ...]. For example, n = 1 
yields the well-known Fibonacci sequence related to the 
golden mean Aa u = [1] = (1 + VE)/2, the case n — 2 cor- 
responds to the octonacci sequence with the silver mean 
A Ag = [2] = 1 + V2, and for n = 3 one obtains the bronze 
mean A Bz = [3] = (3 + y/l3)/22± 

Due to the recursive inflation rule |T|) , these quasiperi- 
odic chains possess a hierarchical structure, which is more 
clearly visible by using the alternative construction rule 
C a = C a -i C a -2 (Cq-i)™ -1 , yielding the same quasiperi- 
odic sequences for a > 2 with Co = w and C\ = s. Fur- 
ther, for given n, the structure of these chains consists of 
only two types of clusters with strong interactions, s n and 
,s' !+1 , which are separated by a single weak bond. The 
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latter property can be related to the eigenstates of the 
system, whereas the hierarchical property has a crucial 
influence on the transport properties for weak coupling. 
We discuss both aspects in detail later. 

To obtain the quantum mechanical eigenstates of these 
systems, we only consider nearest-neighbor hopping and 
hence obtain the tight-binding Hamiltonian 
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H = j2\i)hi+i{i + i\+J2^ £ < 



(2) 



1=0 



1=0 



represented in the orthogonal basis states \l) associated 
to a vertex I. The off-diagonal matrix elements represent 
the kinetic energy in the tight-binding model. The diago- 
nal elements of the Hamiltonian matrix, which represent 
the potential energy of the sites, are set to zero (e; = 0) 
because no energetic disorder is taken into account for 
the studied systems. 

The hopping parameters t are chosen according to the 
letters s (strong) and w (weak) of the quasiperiodic se- 
quence C with t s = 1 and t w = v (0 < v < f). Thus, 
the number of sites of the ath order approximant is given 
by f a + 1. This model can be interpreted as describing 
an electron hopping from one vertex of the quasiperiodic 
chain to a neighboring one, and the aperiodicity is given 
by the underlying quasiperiodic sequence of couplings. 
We then have to solve the discrete time-independent 
Schrodinger equation 

H\^} = E l \¥) 

by diagonalization of the Hamiltonian matrix of Eq. [2j 
Applying free boundary conditions, we obtain f a + l 
eigenstates = Y^i=o 10 ano - the corresponding en- 
ergy values E l . 



III. TIME EVOLUTION OF A WAVE PACKET 
ON QUASIPERIODIC CHAINS 

As outlined in the introduction, the dependency of 
transport properties on the spectral properties of the 
Hamiltonian is not yet fully understood, and thus the 
investigation of transport properties in quasiperiodic sys- 
tems continues to be of special interest. In this section, 
we study the wave-packet dynamics in such systems, in 
particular in the limit of weak coupling where v -C 1. 
The results will then be related to the properties of the 
eigenstates using a perturbation theory approach in the 
following section. 

We investigate the transport properties by means of 
the time evolution of a wave packet |$) = J2i=o®i I0> 
which is initially localized at the center of the quasiperi- 
odic chain, i.e. <!>;(£ = 0) = <5;; with l = \f a /2\. It 
is represented in the basis of the orthonormal eigen- 
states $/(i) = Ei*/o*K*)- The solutions of the time- 
dependent Schrodinger equation then follow by the sep- 
aration approach with ^>\{t) = ^\e~ lE 1 . 
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FIG. 1: Evolution of the width d(t) of a wave packet initially 
localized in the middle of the silver mean chain C^f with 
3364 sites, for several small values of v. The insets show two 
magnified steps for v — 0.03. 



Besides calculating the expansion of the wave packet 
in space, a more detailed analysis can be obtained by 
the computation of the temporal autocorrelation function 
(also known as return probability) 

and the mean square displacement (also called the width) 



d(t) 
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of the wave packet. It is known that a particle's return 
probability decays with a power law C{t) ~ t~ s , where 5 
is equivalent to the scaling exponent of the local density 
of states , 12 ' 13 ! 14 ' 15 and 8=1 refers to ballistic motion. It 
is further known that the spreading of the width d(t) of 
the wave packet shows anomalous diffusion, i.e. d(t) ~ t 13 
with < (3 < i^jj^w^o Here p = o corresponds to 
the absence of diffusion, (3 = 1/2 to classical diffusion 
and (3 = 1 to ballistic spreading. 

In addition, the wave packet dynamics exhibit mul- 
tiscaling, where different moments of the wave packet 
scale with different, nontrivially related exponents 
^21,22,23,24,25,26 \y m i e wav e packet localization implies 
a pure point spectrum, the converse is not true, and 
the more refined notion of semi-uniform localization is 
necessary^ However, the exact relations between parti- 
cle dynamics and singular or absolutely continuous spec- 
tra are less well understood. As a rule of thumb, systems 
with singular continuous spectra exhibit anomalous dif- 
fusion, while absolutely continuous spectra may lead to 
either anomalously diffusive or ballistic dynamics . 21 ' 28 

Consequently, the systems considered here should show 
anomalous diffusion with the mentioned power-law de- 
pendency d(t) ~ t 13 due to their singular continuous 
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FIG. 2: Same as Fig. [T] but for (a) the golden mean model with 4182 sites and (b) the bronze mean model Cg z with 5117 
sites. 
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FIG. 3: Temporal autocorrelation function C(t) of a wave packet initially localized at the center of (a) a golden mean chain 
C\s with 4182 sites, (b) a silver mean chain C^f with 3364 sites, and (c) a bronze mean chain Cf z with 5117 sites. Results are 
shown for small coupling parameters v. 



spectra ; 20 ! 29 which is confirmed by our numerical re- 
sults. Figures [T]and [2] show the development of the mean 
square displacement of the wave packet over time for the 
golden, silver, and bronze mean models in the regime 
of strong quasiperiodic modulation (v < 0.3). While 
the asymptotic behavior can be characterized as anoma- 
lous diffusion, e.g. with scaling exponents not larger than 
/3 Au ~ 0.41, /3 Ag « 0.39, and (3 Bz ~ 0.40 obtained for 
v = 0.2 (for higher values of v and the associated val- 
ues of (3 see Yuan et ali 2 ^) , there are intervals where d{t) 
grows according to a power law d{t) ~ t@ with a constant 
(w-independent) exponent (3 1 which are intercepted by 
flat regimes. As demonstrated by the insets in Fig. [T] in 
these flat regimes d(t) strongly oscillates in a self-similar 
manner, reflecting the hierarchical structure of the sys- 
tem. Nevertheless, the width d(t) remains bounded from 
above by a constant. 

While for the silver mean chain the steps in log d(t) 



have about the same size for a particular value of v, we 
observe small as well as large steps for the golden and 
bronze mean models. This may be caused by the dif- 
ferent underlying inflation rule, which in the case of the 
octonacci chain leads to a symmetric sequence (the se- 
quence is palindromic for this case) and to asymmetric 
ones for the other two systems. However, as a general 
trend we observe an increase of the logarithm of the time 
distance between successive steps with decreasing cou- 
pling constants v in all three models. 

The same behavior can also be observed for the au- 
tocorrelation function C(t), as shown in Fig. [31 All 
three quasiperiodic chains show a stepwise behavior of 
the return probability C(t), where the minimum value of 
C(t) depends on the system size and on the number of 
w bonds. Again the stepwise process is clearly visible, 
where a step here consists of the decrease of C(t) with a 
power law with w-independent exponent 5' , followed by 
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a time interval of constant return probability. The time 
intervals for the power-law behavior and the flat parts of 
both quantities C(t) and d(t) are in correspondence with 
each other. Further, for large time the return probability 
and wave-packet width are bounded by a constant due 
to the finite size f a of the system. 

A more detailed inspection of the wave packet dynam- 
ics reveals that breathing modes are responsible for the 
oscillations, while the wave-packet spreading itself is lim- 
ited to low-amplitude leaking out of the region in which 
it is confined. Eventually, the wave packet expands fast 
to reach the next level of the hierarchy, before the whole 
process repeats. This behavior is even more evident in 
the time evolution of the probability density of such wave 
packets, as shown in Fig.|U At first up to t — e 14 in Fig. [4] 
the wave packet is confined in a narrow range around its 
initial position in the environment of the approximant 
C^ s and oscillates back and forth in this range corre- 
sponding to the strong oscillations of d(t). Then, between 
t = e 16 and t = e 18 the wave packet expands almost bal- 
listically and a significant probability density is found in 
the neighboring C^ s sequences (cp. t > e 18 ). A similar 
behavior can be already seen for the previous levels of 
the hierarchy (cp. the panels up to t = e 10 with those for 
t > e 12 in Fig. [4j) , where the wave packet can be observed 
to spread from the central C^ s structure to a three-fold 
C§ s sequence occurring in the middle of the central C^ s 
chain. 

This indicates that the values of the flat regimes in 
d(t) and C(t) are directly related to the chain structure 
and do not depend on the coupling parameter v. In par- 
ticular, we can give a rough estimate of the correspond- 
ing values of d(t) for the plateaus by assuming that the 
wave packet is uniformly distributed in a confined region 
at the center of the chains. For an approximant C^ s 
the center of the octonacci chain is made up from the 
three-fold sequences C^ g C Ag C^ s with o = a - 2b (b e N, 
< o < a). Confining the wave packet to these re- 
gions we obtain for the approximant C 10 s of Fig. [1] the 
values d(t) = 2.9, 15, 86, 500 for o = 2, 4, 6, 8 respectively. 
Although the wave packet is far from being uniformly 
distributed and some parts of the wave packet can also 
be found outside the confined region (cp. Fig. [4j), our as- 
sumption results in a good reproduction of the plateau 
values observed in Fig. Q] 

This spreading of the wave-packet width to the next 
level can be described by a power law d(t) ~ t 13 with 
exponents f3' An w 0.88, f3' Ag ks 0.85, and j3' Bz » 0.98 de- 
termined for the smallest coupling constant v = 0.03 con- 
sidered here. Computing the results also for all 6 other 
qualitatively different initial positions of the wave packet 
for the octonacci chain we obtained about the same scal- 
ing exponents f3' . For instance for the octonacci chain we 
found P' — 0.81 — 0.85, where some of these differences 
might be caused by fluctuations of the width d(t) which 
are present in the regime of strong expansion and make 
fitting difficult. These scaling exponents tend towards 
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FIG. 4: Snapshots for the evolution of a wave packet initially 
localized at the center of the octonacci chain Cq S with v = 0.1. 
The wave packets are equally scaled for all considered times. 
Additionally, the hierarchical structure of the approximant 
C^ s is visualized by showing the occurrence of the patterns of 
Cg S and C^ s in the chain. 



the exponents obtained for v — > lj2£ which indicates that 
the fast expansion is not governed by the weak coupling, 
but rather a kind of resonance between the different lev- 
els of the hierarchy. Further, for the return probabil- 
ity we find the exponents S' Au w 0.71, S' Ag fa 0.71, and 
S' Bz w 0.76, which are again relatively close to the expo- 
nents for v — > li 20 i 28 However, these values of the scal- 
ing exponent 8' might differ from the exact result, be- 
cause in one dimension we cannot rule out the influence 
of subdominant logarithmic contributions for the consid- 
ered short time intervals in the step-like process^ 

Similar behaviors for d(t) and C(t) have been reported 
before for the Fibonacci chain with strong quasiperiodic 
oscillations ! 28 ! 30 Wilkinson and Austin found the same 
step-like process when studying the spreading of a wave 
packet for Harper's equation of an electron in a mag- 
netic field^ 3 - Based on a qualitative model of the wave- 
packet spreading in the semiclassical approximation and 
on numerical simulations, they argued that a hierarchi- 
cal splitting of the energy spectrum into constant-width 
bands leads to a step- like behavior with /?' = 1, which is 
smoothed due to the (broad) distribution of band widths. 

The value /?' < 1 observed here suggests that there is a 
distribution of band widths in the energy spectrum even 
for v <C 1. The energy levels of the singular continuous 



5 



spectrum for < v < 1 cluster in certain energy ranges 
and form bands with different widths (cp. Sec. lIVp . 
but all the band widths approach zero in the limit v — > 
0i 31 i 32 This suggests that the self-similar spreading of the 
wave packet is only an approximate description of a more 
general multiscale dynamics. 

Further, the self-similarity of quasiperiodic sequences 
was previously used in a renormalization-group perturba- 
tive expansion that provided a great deal of insight into 
the eigenstate properties j 3 - 3 - and showed multiscaling of 
wave packet dynamics. 24 In the following section we fo- 
cus on ergodic rather than hierarchical properties^ By 
an elementary analysis of the perturbation theory of de- 
generate levels at v — for small coupling constants v, 
we show that, on the one hand, eigenstates delocalize for 
any v > 0, in contradistinction to (trivial) localization at 
v = 0. On the other hand, in the limit as v — > 0, eigen- 
states delocalize across only one set of clusters containing 
the same number of atoms, i.e. we obtain a subcluster lo- 
calization due to the breaking of ergodicity. 

IV. RALEIGH-SCHRODINGER THEORY 

Raleigh-Schrodinger theory allows the recursive con- 
struction of matrices in subspaces for a degenerate 
eigenenergy to a given order p, whose diagonalizations 
(the secular problem) yield corrections to the unper- 
turbed eigenenergies. Within this approach, the Hamil- 
tonian is decomposed into an unperturbed system 
and a perturbation i?W with H — + , where 

the hopping parameter v is treated as the perturbation 
yielding = H(v = 0) and = H{v) - H(v = 0). 
Although the accuracy of 0(v p+1 ) of the expansion to 
the pth order is not guaranteed, the theory yields good 
results for small perturbations and preferably large sep- 
arations between the degenerate energy levels. The first 
condition is met since we only consider small values of 
v and the latter one is satisfied since we found that for 
the unperturbed system \E l — E j \ > c with c Au = 0.41, 
c A « = 0.20 and c Bz = 0.11. 

For the considered quasiperiodic sequences the results 
of such perturbation expansions yield two qualitatively 
different types of solutions, depending on the values of p 
and n due to the chain structure mentioned in Sec. [U 
The reason is that the chain consists of strongly coupled 
clusters with n+1 and n+1 atoms, which are weakly con- 
nected via the hopping parameter v. For v = we have 
an unperturbed system with 2n + 3 highly degenerate 
levels, where all eigenstates are localized on individual 
clusters. In higher orders of perturbation theory, these 
localized states then spread first over neighboring clus- 
ters of the same type as the coupling among the clusters 
is taken into consideration, and, for a sufficiently high 
order, delocalize across the whole chain. 

Since the maximal number of letters w between two 
consecutive clusters of length s n+1 is also n+1, the eigen- 
states of these clusters delocalize only in order n+1 of the 



perturbation theory. However, small clusters of length s n 
are connected by at the most 2 (n > 1) or 3 (n = 1) weak 
bonds. More precisely, the dimension of the secular prob- 
lem for each type of cluster separated by not more than 
q weak bonds, changes from 0(q 2 ) at most for p < q to 
0(f a ) for p > q. Only the latter case allows for multifrac- 
tal and/or extended states to be present in the solutions 
of the perturbed system. 

As an illustrating example we analyze the octonacci 
chain with n = 2. In the unperturbed system there are 7 
levels, given by E^l = ±(v% ± l)/2 and = ±y/2, 0. 
Due to the sequence of the ss and sss clusters, the eigen- 
states delocalize only in the 2nd and 3rd order of the 
perturbation theory. In first-order expansion (p = 1), the 
only correction are 6 levels linear in v, splitting off the 
three Eis levels, because -]ssk;ss[- is the only possibil- 
ity for which two clusters of the same type are connected 
by a single w bond. There are also remaining separate 
ss clusters, so that the unperturbed E^ levels are still 
present in the spectrum. For p = 2, all ss clusters be- 
come connected while sss states are still not extended 
due to the existence of -\ssswsswsswsssl- sequences in 
the chain. For p = 3, states belonging to sss levels delo- 
calize as well. Figure [5] shows this splitting of the energy 
levels, where the solutions obtained by different orders 
of the perturbation theory are shown in comparison to 
the energy levels obtained by numerical diagonalization 
of the Hamiltonian H . The results of the perturbation 
approach and the exact energy values are close, although 
only 2nd-order corrections are taken into account. In the 
case of Egs the 2nd-order corrections seem to overcom- 
pensate the error leading to a stronger splitting of the 
energy levels compared to the exact numerical results. 

To investigate the issue of convergence further, we no- 
tice that, even when calculated to all orders of v, the 
secular problems for the two types of clusters give solu- 
tions that are inevitably restricted to the clusters of the 
given type, with zero component on the clusters of the 
other type. For various values of v we check by numerical 
diagonalization whether this can be confirmed. Figure [5] 
shows the total probability that the particle in an eigen- 
state with an energy E l will be on a large cluster 

P s , l+1 (E l ) = J2 (3) 
;es"+ 1 

for the golden, silver, and bronze mean model. Figure 
[6] illustrates how the degenerate eigenvalues for v = 
spread into wider and wider bands with increasing v. The 
results imply that P s n+i strongly depends on the energy 
of the corresponding eigenstate for small v, because it is 
either large for the states belonging to the s n+1 bands 
and vanishes for the states of the s n bands in the limit 
v — ► or shows the reverse behavior. From Fig. [6] we 
can also observe that there are n + 1 bands with high 
probabilities and n bands with low probabilities, in cor- 
respondence with the number of eigenstates generated by 
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FIG. 5: Comparison of energy levels obtained by the perturbation theory approach and by numerical diagonalization of the 
perturbed Hamiltonian H. Results are shown for the octonacci chain with n = 2, a = 6, and (a) ti = 0.1 or (b) v = 0.2. For 
reasons of symmetry only energy values E > are included. 
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FIG. 6: Probability P s n+i(E) that a particle in an eigenstate with a given energy E is on the s n+1 sublattice, for v — 
0.0, 0.1, . . . , 0.9. Only states with E < are shown due to the symmetry of the eigenstate spectrum with respect to E — 0. 
Results are given for (a) the golden mean model Cfg 1 with 6765 sites, (b) the silver mean model with 8120 sites, and (c) 
the bronze mean model Cf z with 5117 sites. Note that there are a few additional states not included in the perturbation theory 
approach, which are caused by the boundary conditions. 



7 



the large and small clusters as obtained by our perturba- 
tion theory approach. However, for v > the probability 
P s n+i is greater than for all energy bands, which im- 
plies that the eigenstates spread over both types of clus- 
ters and thus are ergodic. Here ergodicity denotes the 
spreading of eigenstates over both types of clusters. 

While Fig. [6] shows results for only one iterant of the 
three considered models, we have also compared the re- 
sults for P s n+i of smaller and higher approximants to 
see whether there is any systematic deviation from the 
shown values. We found that, although there are much 
less states for smaller systems or many more states for 
larger approximants, the values of P s n+i do not shift sys- 
tematically for almost all of the states. Instead, the ad- 
ditional points cluster in the same way as in Fig. [6l 

This supports the possibility that eigenstates for an in- 
finite quasicrystalline chain remain ergodic even for very 
small values of v. This is in disgareement with the re- 
sults of the Raleigh-Schrodinger perturbation expansion 
of degenerate levels, which might, nevertheless, still be 
accurate in the limit as v — > since 

. v^a J 1 if E % is caused by s" +1 clusters 

r s n+iyhj ) > < , 

I if E % is caused by s" clusters 

which is a necessary condition but not a sufficient one. 

V. ERGODICITY AND INFLUENCE OF AN 
IMPURITY 

Although the eigenstates remain ergodic even for weak 
coupling v, the wave packet dynamics are strongly influ- 
enced, as outlined in Sec. IIIII This leads to interest- 
ing consequences when a single impurity of strength u 
is placed at a site V by changing a diagonal element of 
the Hamiltonian, i.e. H' = \l')u(l'\. At first we study 
the spreading of wave packets in the presence of such an 
impurity at different types of clusters and then present 
results for the maximum wave packet width when the 
impurity is placed at or near the initial position of the 
wave packet. 

In the first situation we performed several numerical 
experiments for various initial positions of the impurity 
and of the wave packet for different values of u. We found 
that for large u the impurity acts as a barrier, effectively 
cutting the chain into two halves. The consequence is 
that the wave packet is reflected at the impurity indepen- 
dently of its initial site, even if the coupling v is small. 
For u — > 0, on the other hand, the unperturbed wave 
packet propagation of the case u = is restored. 

Understanding the wave packet propagation in the 
regime of intermediate values of u, however, poses sig- 
nificant challenges and surprising results-^- A common 
situation is shown in Fig. [7J for the silver mean model, 
where the wave packet is initially localized at an sss clus- 
ter and the impurity u is placed either on an ss or sss 
cluster near the center of chain. In this case the evolu- 
tion of the wave packet exhibits high sensitivity on the 



position of the impurity, approaching two quite differ- 
ent stationary states. In particular, while in Fig. [7(a) 
the final state is just slightly perturbed from the final 
state for u — 0, in Fig. [7(b) most parts of the wave 
packet are reflected and only a small amplitude can leak 
through the barrier. This kind of wave packet dynamics 
is a consequence of the nearly non-ergodic spreading of 
the eigenstates for the two different types of clusters for 
small v as discussed above (cp. Fig. [6]). The explanation 
is that the wave packet is constructed by a superposition 
of all eigenstates and because it is initially localized on 
an cluster, eigenstates caused by this type of cluster 
are contributing with a much higher probability than the 
eigenstates of the s n clusters. Consequently, the impu- 
rity is felt as a barrier by the wave packet when placed 
on the same type of cluster as the initial position of the 
wave packet. 

In the second situation we address the influence of an 
impurity placed at or near the initial site of the wave 
packet on the dynamics by studying the dependence of 
the final wave packet width on the impurity strength u. 
Figure [5] shows the maximum value of d(t) attained in 
the course of the evolution of a wave packet in the oc- 
tonacci chain, which was initially localized at the site 
'o = r/a/2] ■ Here only systems C^ 6 with odd a and with 
the impurity placed at the initial site of the wave packet 
or its left neighbor site are shown. To obtain these max- 
imum values of the width we perform the calculations as 
in Fig.[TJ but now for the perturbed system and different 
values of the impurity strength u up to very large times, 
where the system is governed by finite size effects and 
d(t) becomes constant. 

The results show that, for small u, the width of the 
wave packet equals the results for the unperturbed sys- 
tem for both positions of the impurity. For large u, we 
obtain a strongly localized final wave packet when the 
impurity is added at the initial site Iq of the wave packet 
and a constant width of the wave packet when placing 
the impurity at the left neighbor site Iq — 1. In the first 
case, the expansion of the wave packet in the eigenstate 
basis is dominated by strongly localized wave functions 
caused by the large impurity u and the wave packet can 
no longer spread across the chain. In the latter case, 
the impurity acts as a barrier placed at the center of the 
chain, and consequently the wave packet is reflected and 
only spreads across one half of the system as in Fig. [7(b) , 
and d(t) reaches a plateau in Fig. [8] Nevertheless, the 
width d(t) is reduced compared to an unperturbed sys- 
tem with half the system size because in the presence 
of an impurity always some localized eigenstates, which 
do not spread across the quasiperiodic chain, occur and 
contribute to the expansion of the wave packet. 

However, in between these two extremes there is a wide 
range of values of u for which the final width of the wave 
packet is significantly reduced even for u <C v, signal- 
ing dynamical localization. There are nevertheless sev- 
eral well-defined peaks in Fig. [5] for some values of u at 
which the maximum wave-packet width is significantly 
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FIG. 7: Snapshots of the evolution of two wave packets in the presence of an impurity for the octonacci chain with n — 2, a — 8, 
and v = u = 0.1. The wave packet is initially localized on a large cluster in a local environment -]ws x ssw[- at the site x. In the 
two panels (a) and (b) the impurity is located on a small cluster or a large cluster in a local environment ^wSbSSWSaSwssswl- 
at the sites indicated by a and 6, respectively, as visualized by a vertical line in each panel. For easier comparison, the vertical 
dashed line in each panel marks the position of the impurity in the other panel. The long-time wave packet dynamics exhibit 
high sensitivity on whether the impurity is located on the same type of cluster or not as the one on which the wave packet was 
initially localized. 
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FIG. 8: Maximum width of a wave packet attained during its 
evolution in the presence of a single impurity of strength u 
for the octonacci chain with n = 2 and v = 0.1. The wave 
packet is initially localized at the center of the chain x = lo in 
the local environment ■]wssws ui s x swssw[- and the impurity 
is either placed at the same site mo = x or its left neighbor 
site Ui. Note that we only included approximants C^ g with 
odd a because for even a the local environment of the initial 
position of the wave packet is different. 



enhanced, compared to the cases of slightly smaller and 
slightly larger values of u. These peaks persist for differ- 
ent system sizes and for different positions of the impurity 
U, although the positions and structure of the peaks can 



change. Figure [5] shows that the peaks for the four sys- 
tems considered there occur at the same strength of the 
impurity, especially for relatively large values of u. 

Having a closer look at the eigenstates caused by the 
impurity, we found that those peaks appear at impu- 
rities u, where at least some of these perturbed states 
coincide with the bands in the energy spectrum of the 
quasiperiodic approximants shown in Fig. [SJ In this case 
the states caused by the impurity hybridize with the un- 
perturbed states of this band and consequently the wave 
packet is able to spread along the chain easier. This also 
explains the differences of the peaks for a = 7 and a = 9 
in Fig. [8j because the structure of the energy spectrum 
of the 9th approximant is more complicated and thus 
additional peaks occur. 

Further, in Fig. [5] we compare the maximum widths 
of the wave packet for impurities placed at the 1st left 
neighbor and the 3rd left neighbor in systems C^ s with 
even a. It is clearly visible that the maximum widths of 
the wave packet for these systems are almost identical for 
the two different positions of the impurity. The reason 
is that the impurities are in both cases located at the 
edges of the ss cluster, which is adjacent to the cluster, 
where the wave packet is initially localized. In contrast 
the curve of the approximant Cg g shows a completely dif- 
ferent behavior, because here the wave packet is located 
in a different local environment in the beginning. 

Further, by repeating similar numerical experiments 
for various values of the coupling strength v, we found 
that the peak structure becomes less distinct with in- 
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FIG. 9: Same as Fig. [8] but for even a with the impurity of 
strength u located at the 1st and 3rd left neighbor. The wave 
packet is initially localized at the center x = lo of the chain 
in the local environment ^wsssWu^sSmWxSSWSSSwl- . 

creasing v. For the silver mean model we observed that 
it persists up to v » 0.4. For this value of v the widths 
of the energy bands become smaller than the gaps be- 
tween them, which in turn means that the impurity- 
related eigenstates coincide very often with the energy 
bands and thus peaks merge so that almost no valleys 
occur. 

These results show that in quasiperiodic quantum 
wires one can strongly influence the long-range electronic 
transport properties by inducing local perturbations at 
different positions and of various strengths. The char- 
acteristics are related to the nature of the eigenstates, 
which spread only across one type of cluster in the limit 
v — > 0. Knowing the structure of the energy bands and of 
the eigenstates allows one to design quasiperiodic chains 
with impurities that can act as sort of control gates. 



with time intervals of power-law growth, followed by a 
regime with confined wave-packet width. Nevertheless, 
the average wave packet dynamics can be classified as 
anomalous diffusion. These results are consistent with 
the literaturej 20 ! 23 ' 28 The stepwise behavior is caused by 
the hierarchical structure of the chains, leading to the 
confinement of the wave packet until it expands fast to 
reach the next level of the hierarchy. 

The perturbation theory approach allowed us to draw 
a connection between the structure of the weakly coupled 
clusters and the localization characteristics of the eigen- 
functions 'J, which only become delocalized in sufficiently 
high orders of the perturbation expansion. This happens 
when clusters of a specific type become connected, i.e. in 
order n + 1 of the expansion for large clusters and in the 
2nd or 3rd order for small clusters. 

However, while we obtained that the eigenfunctions 
spread ergodically over all clusters of the chain for v > 0, 
in the limit in we found that ergodicity is broken 
and the eigenfunctions can only spread across clusters of 
the same type. The latter case has a significant influence 
on the long-term wave packet dynamics in the presence 
of local perturbations. We showed that the initial site 
determines whether the wave packet is reflected at the 
impurity or whether it can tunnel through the impurity, 
and that the maximum width of the wave packet is af- 
fected by the strength and the position of the impurity. 
This behavior might be used to construct control gates 
for the electronic transport in quasiperiodically modu- 
lated quantum wires. 

Further research is planned in the following directions: 
It is intended to consider the properties of more complex 
systems, as e.g. the associated labyrinth tilings in two 
and three dimensions. Further, instead of investigating 
an on-site disorder it would be also interesting to study 
the influence of phason disorder caused by exchanging 
certain strong and weak bonds. 



VI. CONCLUSION 

In this paper, we considered the electronic transport 
in one-dimensional quasiperiodic chains consisting of two 
types of clusters which are weakly coupled by a hop- 
ping parameter v. The investigations of the wave-packet 
dynamics revealed the occurrence of a stepwise process 
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